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NLO perturbative QCD predictions for 77 -> M + M~ (M = it, K) 
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We report the first complete leading-twist next-to-leading-order perturbative QCD predictions for 
the two-photon exclusive channels 77 — > M + M~ (M = ir, K) at large momentum transfer. The 
asymptotic distribution amplitude is utilized as a candidate form for the nonperturbative dynamical 
input. Comparison of the obtained results with the existing experimental data does not provide 
sufficiently clear evidence to support the applicability of the hard-scattering approach at currently 
accessible energies. 
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Exclusive processes at large momentum transfer pro- 
vide a particularly suitable testing ground for pertur- 
bative QCD (PQCD). In the hard-scattering approach 
(HSA) 1] the amplitudes of these processes are given 
in terms of a convolution of a hard-scattering amplitude 
of collinear hadron constituents, and distribution ampli- 
tudes summarizing soft physics. Although the HSA is 
widely used, its applicability at currently accessible en- 
ergies is still being questioned. 

The leading-order (LO) PQCD predictions have been 
obtained for many exclusive processes, but owing to their 
sensitivity to the choice of the renormalization scale and 
scheme, they do not have much predictive power. To 
stabilize the LO results, and to achieve a complete con- 
frontation between theoretical predictions and experi- 
mental data, it is crucial to obtain higher-order contri- 
butions. However, owing to great technical difficulties 
involved, the complete next-to-leading-order (NLO) pre- 
dictions were obtained only for the pion electromagnetic 
form factor 0, Q and the pion transition form factor 
ill 

The two-photon annihilation into flavor-nonsinglet 
helicity-zero mesons, 77 — » M + M~ (M = n, K), is one 
of the simplest hadronic processes. Owing to the point- 
like structure of the photon, the initial state is simple, 
with strong interactions present only in the final state. 
The energy behavior, angular behavior, and normaliza- 
tion of the cross section of this process are nontrivial 
predictions of PQCD. This process has been the subject 
of a number of experimental investigations. The most 
recent data (DELPHI 0, ALEPH @, and Belle 0) are 
obtained for the center-of-mass energies W < 6 GeV and 
the scattering angles in the range | cos6>| < 0.6. 

In the HSA, at large momentum transfer (high W and 
large 9), the process is described by the helicitiy ampli- 
tude 
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xT H (\\';x,y;W,6;n 2 F ,fi 2 R ), (1) 



The function $m{%, Mf) is the meson distribution am- 
plitude (DA), representing the amplitude for the meson 
consisting of a qq pair, with the quark and antiquark 
collinear and on shell, and sharing fractions x and 1 — x 
of the meson's total momentum. There is a long on- 
going debate in the literature Q regarding the correct 
form of the DA at currently accessible energies. For 
the purpose of this work we approximate the DA by 
its asymptotic form and choose §m{x, ix 2 p) — ^%f(x) = 
/m\/3/2 x(l — a:), where /m is the meson decay constant 
(U(K) =0.131 (0.160) GeV). 

The function T H (XX';x,y; W, 9; fip, fi 2 R ) in Eq. (JTJ is 
the hard-scattering amplitude for producing collinear me- 
son constituents from the initial photon pair, is free of 
collinear singularities, and has a well-defined perturba- 
tive expansion. 

At the NLO the helicity amplitude in Eq. (JTJ is given 
by the NLO result for Th, obtainable from the LO and 
the NLO contributions to the 77 — > (qq) + (qq) amplitude, 
with massless quarks collinear with outgoing mesons. 
These contributions arise from 20 tree and 422 one-loop 
Feynman diagrams, respectively. Distinct diagrams are 
shown in Fig. ^ The other diagrams are obtained from 
these using various symmetries. 

The LO prediction for the process 77 — ► M + M~ (M — 
7T, K) is due to Brodsky and Lepage (BL) Making 
use of the similarity of the LO expressions of the helicity 
amplitudes for the process in question and the timelike 
meson electromagnetic form factor Fm, they obtained 
the approximate relation 
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where A and A' are photon helicities, /i F is the factoriza- 
tion and the renormalization scale. 



Despite the fact that the complete LO prediction is given 
in Ref. , the approximate relation of Eq. J3J) is what 
is usually referred as the Brodsky-Lepage prediction. 

The first attempt to perform the NLO analysis was 
done by Nizic |l0j. Lacking powerful analytical meth- 
ods for calculating one-loop Feynman diagrams with five 
and six internal lines, the NLO prediction in Ref. |Hj 
was obtained with the simplest equipartition meson DA 
<&m oc d (x — 1/2). Also, the contribution of the diagrams 
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LO diagrams: 



f2^\ D 



NLO diagrams: 
group A group B 



Al^-- A2^\ A3^-- A4^^ A5 



B3 B4 B5 



A6^-^_ A7^^ A8^-~^ A9^^ AltJ^^ B(i B7^^^ B8^-^ B9^--, BIO 



All^^ A12~~-~, AliT-^ Al4^^ A15 



Bll^^ B12^^ B13^^ B14\, B15 



group C 



group D 



Cl""---. C2^-^ C3^-^ C4^--. CS 



Dl~\. D2^--. D3^-^ D4^^ D5 



Dti""---. D7^-~, D8^\ D9^^ DltT^ 



group E 



El^^ E2 



FIG. 1: Distinct diagrams contributing to the 77 — > (qq) + 
(qq) amplitude at the NLO. The total number of LO diagrams 
is 20, and the total number of NLO diagrams in groups A, B, 
C, D, and E are 96, 156, 94, 70, and 6, respectively. 



of g roup E was omitted. Therefore, the prediction of Ref. 
cannot be considered complete. 

Here we present the results of a complete NLO anal- 
ysis of the process 77 — > M + M~ (M — tt, K). Using 
the method of Refs. [llj, we evaluated all the one- loop 
diagrams analytically. It should be pointed out that the 
computation of one- loop diagrams with five and six inter- 
nal lines represents a highly demanding task. The details 
will be given elsewhere. 

Our complete NLO PQCD result for the cross section, 
in the angular range | cos#| < 0.6, assuming the asymp- 
totic DA and using the MS renormalization scheme, is 



cr(77 -*-M + M") = 



4 1-035 
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W 2 



,(3) 



where a s (n 2 R ) = 47r/( /3n ln(/^/A 2 ) ) is the coupling con- 
stant renormalized at ^ 2 R (function of W 2 only), with 
Pa = 11 - 2n//3 (n f = 3) and A = 0.2 GeV. As a con- 
sequence of using the asymptotic DA, the result © docs 
not depend on ^l 2 f . The corresponding result for the dif- 
ferential cross section is too lengthy to be presented here. 

As it is seen from Eq. J2J, the NLO corrections, at 
accessible energies (W < 6 GeV), are substantial. Thus, 
for = W 2 (pragmatical choice), the ratio of the NLO 
to the LO contributions is 0.9 for W = 4 GeV, where 
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FIG. 2: Cross sections for the processes 77 — > M + M 



a s (/i|j) — 0.23, and becomes less than 0.5 only for W > 
45 GeV, where a s (fj, R ) = 0.13. Obviously, the problem 
of convergence relates to the size of the NLO coefficient, 
not to the coupling constant. 

It is important to observe that the large contribution 
to the NLO coefficient in Eq. comes from the /?o term, 
arising from the vacuum-polarization diagrams, and it is 
desirable to resum them into the running coupling con- 
stant. This can be achieved using the Brodsky-Lepage- 
Mackenzie (BLM) scale setting procedure |l^ . 

However, the fact that the form of the effective cou- 
pling at low energies, as well as its behavior at timelike 
scales are not determined, together with the fact that 
the process at hand contains both spacelike and timelike 
BLM scales, does not allow one to perform the BLM pro- 
cedure in a completely satisfactory way. For the sake of 
clarity, we discuss only the BLM improved LO predic- 
tion LO(BLM), which is of the broadest interest. The 
BLM scales were determined for each LO diagram sep- 
arately, for different helicities, angles, and quark mo- 
menta. The mean- value theorem for integration was used 
to avoid too low scales in the integration indicated in Eq. 
Q). To deal with the timelike scales, the continuation 
ln(/z 2 /A 2 ) — > ln(/i 2 /A 2 ) — iir was used in the coupling 
constant. 

In Figs. [21 and [21 we compare our predictions with 
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FIG. 3: Angular dependence a 1 d<r/d|cos#| for the 77 — > 
M + M~ processes. The experimental points are from Ref. Q- 



the recent experimental data and the BL estimate in 
Eq. J2J). To assess the theoretical uncertainty related 
to the choice of the scale ^ 2 R in Eq. J3J), in Fig. |3 
we show the range of the NLO predictions (shaded re- 
gion) obtained by varying [i 2 R from fi R = W 2 down 
to the scale where the NLO prediction reaches maxi- 
mum, /i R « W 2 /lb (the scale determined by the fastest- 
apparent-convergence principle [13|). The NLO predic- 
tion (not shown in Fig. |5J for fi R W 2 /20, the scale 
set by the principle of minimal sensitivity 14], is lo- 
cated slightly below the upper limit. The solid line in 
Fig. [21 is our LO(BLM) prediction, while the dotted 
line is the BL estimate based on Eq. Q with the form 
factors approximated by F n (W 2 ) » 0.4 GeV 2 /W 2 and 



F K {W 2 ) « f%lBFv{W 2 ) = 0.6GeV7^ 2 , as in the 
original paper J9J. However, the recent measurements 
EH, suggesting that F^ [K) {W 2 ) « 1.01 (0.85) GeV 2 /W 2 
(around W 2 = 13.48 GeV 2 ) put the success of the widely 
cited BL prediction (dotted line) to question. Figure [5] 
shows that our results for the cross section are still sig- 
nificantly, but less than an order of magnitude, below the 
data. 

Next consider the W dependence of the cross section. 
Using the parametrization (7(77— >M + M~) oc W~ nM the 
results obtained by Belle give n T = 7.9 ± 0.4 ± 1.5 and 
n K = 7.3±0.3±1.5,for 3.0 GeV < W < 4.1 GeV. On the 
other hand, we find from Eq. © that the power n„(= 
n K ) takes the values 6.9(7.4) for n\ = W 2 {W 2 /lb). 
The LO(BLM) and BL predictions are 6.7 and 6, re- 
spectively. As for the ratio of kaon to pion cross sec- 
tions, our prediction equals f^l ft = 2.23, coinciding 
with the BL prediction, while the Belle results for 
3.0 GeV < W < 4.1 GeV give 0.89 ± 0.04 ± 0.15. 

In Fig. |21 we show the angular dependence of the ratio 
a- 1 da/d\ cos0|. The range of the NLO and LO(BLM) 
predictions is obtained by varying W from 2.4 GeV to 4.1 
GeV (range covered by the Belle experiment 0). The BL 
curve, based on Eq. J2Jl, corresponds to 1.227/ sin 4 8. To 



keep figure synoptic, experimental points are given only 
for W = 3.2- 3.3 GeV (middle of the Belle range 0). 

As it is seen from Figs. El and the energy and the 
angular dependence of our NLO and LO(BLM) predic- 
tions for the cross section are in very good agreement 
with the data. As for the absolute normalization, our 
results are still significantly below the data. Obviously, 
the normalization represents a problem. 

This problem has been the subject of a long debate, 
and critics use it as an argument to discard, at exper- 
imentally available energies, either the HSA or the 
asymptotic DA as an appropriate distribution amplitude 
|17j . The strength of criticism depends on the size of 
higher-order corrections. Here we have shown that the 
NLO corrections can be substantial, so a part of the nor- 
malization problem can be attributed to the size of un- 
calculated higher-order corrections. However, assuming 
reasonable scale variation of the NLO prediction, or con- 
trolled perturbative series (LO > NLO term > NNLO 
term) it is not to be expected that inclusion of higher- 
order corrections would lead to a solution of the normal- 
ization problem. On the other hand, good agreement 
between theory and experiment for angular dependence 
(Fig. |3J, and the disagreement for the ratio of the kaon 
to pion cross sections suggest that the specific form of 
the DA certainly plays a relevant role. An analysis of 
various choices of the DA, however, is outside the scope 
of this letter. 

As a quantity which is expected to be largely insen- 
sitive to both the form of the effective coupling as well 
as to the choice of the DA, and as such representing a 
nontrivial PQCD prediction, we next consider the ratio 
of (7(77 -> M+M-) and \F M {W 2 )\ 2 at the NLO. 

The result of our NLO calculation of the timelike me- 
son form factor is of the form 



I^m(^ 2 )| 2 = /. 
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and is consistent with the NLO result for the spacelike 
form factor |2|. On the basis of Eqs. © and J3J, and 
using the same we find 



W 



, (7(77^ M + M~) 
\F M (W 2 )\ 2 



= 638 [1 - a s (/4) 0.306] nbGeV 2 



(5) 

A comparison of this result with the only two reliable 
experimental points, based on the data from and |l5| . 
is given in Fig. 0] The solid line corresponds to Eq. JSJ, 
with fi R = W 2 , the dotted line is the BL estimate fol- 
lowing from Eq. the dot-dashed line is the complete 
LO prediction (equal to 638 nb GeV 2 , according to Eq. 
©), while the dashed line is the ratio of the LO(BLM) 
results with the BLM procedure performed separately for 
the cross section and the form factor. As it is seen, our 
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FIG. 4: Ratio W 2 (7(77 -> M+M - VI-Fm {W 2 )?. The exper- 
imental points are obtained from the data of and [TEjl . 



results for the ratio are above the data by approximately 
a factor of 2 to 3 which is in much better agreement with 
the data than the prediction for the cross sections based 
on Eq. (PJ). As Fig. 0]shows, the usage of the BLM scales 
leads to much better agreement with the data. This, to- 
gether with the fact that the LO(BLM) prediction lies in- 
side the range of the NLO prediction (shaded area in Fig. 
EJ), supports the idea behind the BLM procedure accord- 
ing to which the BLM determined scales represent the 
relevant physical scales (putting aside understatements 
regarding the effective coupling) . 

In conclusion, the NLO results obtained for the process 
77 — * M + M~ (M = it, K) are not conclusive enough 
to allow a definite statement regarding the applicability 
of the HSA at the currently accessible energies above 3 
GeV. However, improvements are possible in two direc- 
tions: the first, the more appropriate choice of the DA 
(possibly in the direction suggested in Ref . ^3] ) , and the 
second, a better understanding of the resummation for 
the processes containng both spacelike and timelike scales 
(possibly in the directions suggested in Refs. (lilH^ |). 

The authors would like to thank S. J. Brodsky and 
N. Kivcl for useful comments. 
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